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' . 1 ' Abstract 

\Q ' Certain aspects of the behaviour of the gravitational field near null and spatial infinity for 

, the developments of asymptotically Euclidean, conformally flat initial data sets are analysed. 

Ideas and results from two different approaches are combined: on the one hand the null infin- 
ity formalism related to the asymptotic characteristic initial value problem and on the other 
the regular Cauchy initial value problem at spatial infinity which uses Friedrich's represen- 
tation of spatial infinity as a cylinder. The decay of the Weyl tensor for the developments of 
the class of initial data under consideration is analysed under some existence and regularity 
assumptions for the asymptotic expansions obtained using the cylinder at spatial infinity. 
Conditions on the initial data to obtain developments satisfying the Peeling Behaviour are 
5^ , identified. Further, the decay of the asymptotic shear on null infinity is also examined as one 

approaches spatial infinity. This decay is related to the possibility of selecting the Poincare 
group out of the BMS group in a canonical fashion. It is found that for the class of initial 
• i-h , data under consideration, if the development peels, then the asymptotic shear goes to zero 

■ at spatial infinity. Expansions of the Bondi mass are also examined. Finally, the Newman- 

Penrose constants of the spacetime are written in terms of initial data quantities and it is 
shown that the constants defined at future null infinity are equal to those at past null infinity. 

PACS: 04.20.Ha, 04.20.Ex, 04.20. Gz 



o 
o 



1 Introduction 

Penrose's definition of asymptotic simplicity is an attempt to provide a characterisation of isolated 
systems in General Relativity [571 [551 [55] , It offers a framework on which diverse notions of 
physical interest can be defined and handled in a precise and rigorous manner. There is, however, 
a notorious lack of nontrivial examples where the associated formalism — see e.g. [26, 36, 39ll55j— 
has been actually been used to extract some physics. The obvious reason for this is the scarcity of 
radiative exact solutions to the Einstein field equations other than the boost-rotation symmetric 
ones [2J. A more conspicuous consideration in this regard is the fact that radiative spacetimes 
have to be constructed starting from an initial value problem where, say, some Cauchy initial 
data are provided. Up to fairly recently, there has not been a way of relating the properties of 
the initial data to the radiative properties of the development. For example, how does an initial 
data set have to be so that the resulting spacetime peels'! — that is, the components of the Weyl 
tensor in an adapted gauge having a distinctive decay in terms of powers of an affine parameter 
of the generators of outgoing light cones [351 H2] • Another natural question on these lines would 
be how does the ADM mass and the Bondi mass relate? Similarly, is there any relation between 
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the angular momentum denned at spatial infinity and that defined at null infinity? Which classes 
of initial data allow to select the Poincare group out of the group of asymptotic symmetries — the 
BMS group — in a canonical fashion? 

Arguably, there is a shortage of general results about the evolution of initial data sets for 
the Einstein field equations near spatial infinity. This particular has made the idea of relating 
the structure of initial data with radiative properties a daunting endeavour. However, work by 
Friedrich 20, 23] on the regular initial value problem at spatial infinity for time symmetric initial 
data sets has provided a tool to address questions similar to the ones raised in the first paragraph 
through a systematic approach. Friedrich has introduced a representation of spatial infinity as a 
cylinder — the cylinder at spatial infinity, X — in stark contrast to the usual representation as a 
point. The cylinder at spatial infinity can be regarded as a limit set of incoming and outgoing 
characteristics of the Einstein field equations, and as such it happens to be a total characteristic. 
As such, it allows to transport information from the Cauchy surface to null infinity without being 
contaminated by any sort of boundary conditions. 

The construction which led to the cylinder at spatial infinity allows an unfolding of the evo- 
lution process, which in turn can be analysed in all detail and to the desired order. The latter 
can be thought of as enabling the construction of a certain kind of asymptotic expansions which 
are completely determined by the Cauchy data near infinity. Friedrich's seminal work has been 
extended in several directions, automatising, in some sense, the calculation of the asymptotic 
expansions of the relevant field quantities (50] [49] and to include more general classes of initial 
data sets (52] [53]. The analysis described in these references has exhibited the existence of a 
hierarchy of obstructions to the smoothness of null infinity which goes beyond that arising from 
a mere consideration of the linear aspects of the field equations as described in e.g. [20] [46]. The 
existence of the hierarchy of obstructions has nurtured a rigidity conjecture on the smoothness 
of null infinity which, broadly speaking, states that the only Cauchy initial data sets which have 
a development with a smooth null infinitjQ are those which are stationary close to infinity. Sim- 
ilarly, a given degree of differentiability at null infinity, say C k , would only be possible if the 
initial data is stationary to a given order — which depends on k. The kind of initial data required 
can be constructed by means of the gluing techniques of (9] [TO] [7] . Intuitively, the absence of 
gravitational radiation near spatial infinity which the rigidity conjecture asserts, suggests that the 
incoming radiation has to die off in the infinite future, and at the same time the system cannot 
have emitted gravitational waves for all times in the infinite past. 

Now, the calculation of the asymptotic expansions which arise from the construction of the 
cylinder at infinity is performed in a gauged which is adapted to a Cauchy initial value problem 
— remarkably, the gauge also allows to read off the structure and location of infinity directly from 
the initial data. On the other hand, the discussion of the gravitational field near null infinity is 
usually done using a gauge which is based on geometric structures of null hypersurfaces. Having 
two different gauges, hampers, to a certain extent, extracting the physics of the system. It also 
difficults the attempts to assess the relevance of the presence of obstructions to the smoothness 
to null infinity in the asymptotic expansions. One of the aspects of the analysis of the asymptotic 
behaviour of the gravitational field via the definition of asymptotic simplicity and related con- 
structions is that it allows to rephrase the decay of the fields in terms of differentiability at the 
conformal boundary. Thus, it is natural to expect that the presence of obstructions to smoothness 
in the asymptotic expansions may result in a modified peeling behaviour — this idea has already 
been explored in [51j . 

A first study of the transformation between the gauge used in Friedrich's analysis of spatial 
infinity and the gauge generally used to discuss null infinity has been given in [24] . In the present 
article their analysis will be extended to the developments of a class of initial data sets which 
are not time symmetric. The consideration of initial data sets with a non-vanishing second 
fundamental form will allow us to examine some time asymmetric aspects of the structure of 
the conformal boundary of asymptotically flat spacetimes which up to now have remained terra 
incognita. On the other hand, the discussion shall be restricted to conformally flat initial data 

iflere and in the sequel, smooth will always mean C°°. 

2 In this article, a gauge choice is understood as a certain choice of coordinates, vector or spinorial frame, and 
also a choice of representative in the conformal class of the physical spacetime. 
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sets. This is done for conciseness and for the ease of calculations and presentation. Previous 
analysis with conformally flat initial data sets — see e.g. [501 H2] — have shown that most of the 
crucial phenomena and structure observed in this setting admit a counterpart when discussing 
non-conformally flat data sets — like in [49] [53]. Furthermore, the detailed understanding of 
conformally flat initial data sets is of relevance in view of their use in the numerical simulation 
of black hole spacetimes — see e.g. [27] for a recent discussion on this. A cautionary note is, 
however, due: most of the conformally flat initial sets considered in the numerical simulations are 
boosted. This is done on physical grounds. Here, initial data sets with linear momentum — i.e. 
boosted — will not be considered. The reason being that, as it will be discussed in the sequel, 
they are more technically involved and less smooth. The effects of this lower regularity are of 
interest in themselves and will be discussed elsewhere. 

The article is structured as follows: section [2] discusses some ideas of the basis set up and fixes 
some notation; section [3] describes some aspects of the "cylinder at spatial infinity" -formalism; 
section [J] considers relevant aspects of the class of initial data sets to be used; section [5] muses over 
the asymptotic expansions that can be calculated using the cylinder at spatial infinity; section [6] 
goes briefly about the construction of a frame adapted to null infinity — the NP gauge — while 
section [7] discusses how it is related to the gauge used in the expansions of section [5] This 
article revisits the analysis carried out in [24] of "Bondi type" systems — what here is called the 
Newman-Penrose gauge — and combines it with the asymptotic expansions for the development 
of conformally flat initial data sets which have been calculated in [52] . In section [8] analysis 
of the implications of these expansions and the related obstructions to the smoothness of null 
infinity is done — under some existence and regularity assumptions. In section [9] the machinery 
is used to calculate expansions of the spin coefficient a on null infinity in the NP gauge. The 
behaviour of this coefficient is related to various physical aspects of the theory of isolated systems. 
The implications of these expansions are analysed. It is shown that if the spacetime peels, then 
the behaviour of a is such that it is always possible to single out the Poincare group out of 
the asymptotic symmetry group in a canonical way. Expansions of the Bondi mass are also 
considered. In section [10] the formalism is used to write down the Newman-Penrose constants of 
the developments of the class of initial data under consideration in terms of initial data quantities. 
These expressions allow to show that the constants at future null infinity are equal to those at 
past null infinity. 

There is much more to the aforementioned analysis than what is presented in this article. 
In particular, most of the calculations involved in the analysis have been carried out using ex 
profeso scripts in the computer algebra system Maple V. A first subset of the scripts constructs 
the initial data; another subset determines the coefficients of the asymptotic expansions in terms 
of initial data quantities; finally a third subset calculates the transformation between gauges and 
allows to calculate the asymptotic expansions of NP objects. In order to keep the presentation 
at a reasonable length, the explicit expressions of intermediate steps will be omitted, and more 
emphasis will be put on to the qualitative aspects of the results. Some familiarity with the use 
of spinors — in particular the space-spinor formalism of [41] — will be assumed. Details of the 
diverse other formalisms used — the cylinder at spatial infinity, the formalism of null infinity, the 
construction of gauges, and how to relate them — is kept to the minimum necessary to motivate 
the analysis and to point out the many subtleties arising. In any case, the reader is referred to 
the various references cited for more complete details. Finally, are some concluding remarks and 
a brief appendix containing some useful spinorial expressions. 

2 Basic setup 

Let (.M,<? M „) be a vacuum spacetime arising as the development of asymptotically Euclidean 
Cauchy initial data (S,h a /3,Xa0)- Let C >1 uXp be the conformally invariant Weyl tensor of the 
metric g^ v . Tilded quantities will refer to quantities in the physical spacetime, whereas untildcd 
ones will denote, generically, quantities in an unphysical — i.e. conformally rescaled — spacetime. 
Indices //, z/, A, . . . (second half of the Greek alphabet) are spacetime indices taking the values 
0, . . . , 3; the indices a, f3, 7, . . . are spatial ones with range 1, 2, 3. The Latin indices a, 6, c, . . . 
will be used in spatial expressions which are valid for a particular coordinate system (usually a 
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Cartesian normal one). The indices k, . . . are spatial frame indices, while A, B, C, .. . will be 
spinorial indices taking the values 0, 1. 

The initial data sets are taken to be maximal — that is, \ = X° Q — 0, and the metric to be 
h a /3 negative definite — the latter in order to ease the later use of spinors. For simplicity and 
dcfinitcness, it shall be assumed that the initial hypersurface contains only one asymptotic end 
— that is, a region which is diffeomorphic to R 3 minus a ball. 

The analysis shall be concentrated in the region of spacetime near null and spatial infinity. 
In order to discuss the behaviour of the development of the initial data sets in the asymptotic 
region use will be made of the so-called conformal picture. Let i denote the point at infinity 
corresponding to the asymptotic end under consideration. The point i is obtained by conformally 
compactifying the hypersurface S with a conformal factor f2 which can be obtained by solving 
the Einstein constraint equations — see section 0$. Let S be the resulting compact 3-dimensional 
manifold. The manifold is identified in a standard way with <S\ {i}. Most of the discussion will be 
carried out in a sufficient small neighbourhood B a (i) of radius a centred on i and on the resulting 
associated domains of influence and dependence Af = J + {B a {i)) U J~(B a (i)). 

As it is customary, the discussion of the initial data will be carried out in terms of the tensors 

h a (3 = Q 2 h a p = ■&~ 4 h al 3, Xa(3 = = ^~ 2 Xa/3, 4>af3 = ^~ 1 Xa/3 = ^Xap (1) 

on the conformally rescaled manifold S. In this work attention will be restricted to conformally 
flat initial data sets and their developments. This assumption will prove fruitful and will allow 
to infer features of the development of more general types of initial data sets. Under these 
assumptions coordinates x a can be chosen such that: 

h = -5 ab Ax a &x\ in B a {i). (2) 

Unless otherwise stated, {x a } denote some normal coordinates with origin at i based on a h- 
orthonormal frame ej. 

3 The Friedrich-gauge 

In [50] a certain representation of the region of spacetime N = J + (B a (i)) U J~(B a (i)) — which 
can be thought of as that close to null infinity, J 1 = U and spatial infinity, i° — has been 
introduced. The standard representation of this region of spacetime depicts i° as a point. In 
contrast, the representation of |20j depicts spatial infinity as a cylinder — the cylinder at spatial 
infinity. The technical and practical grounds for introducing this description have been discussed 
at length in the seminal reference. The original construction in [5U] was carried out for the class 
of time symmetric metrics with analytic conformal metric h a b. However, as seen in [52[ 153] . the 
construction can be adapted to settings without a vanishing second fundamental form. 

Friedrich's construction makes use of a blow up of the infinity of the initial hypersurface 
to S 2 . The blow up of i to S 2 requires the introduction of a particular bundle of spin- frames 
over B a (i). Let taa' be the spinorial version of a timelikc vector r M . Assume the normalisation 
t~aa'T AA = 2. Let SU(B a (i)) denote the fibre bundle of spinorial dyads {Sa}a=o,i — {0,4,6,4} 
over B a {i) with projection 71", such that taa' = oaoa 1 + i>A^A J ■ The fibres of SU{B a (i)) are of type 
5f/(2,C) = {t A B e SL(2,C) I T A A't A B t A ' B , = t BB '}- With this notation, define the set 

C a = {6( P ,t A B ) € SU{B a {i)) I |p| < a, t A B e SU(2,C)}, (3) 

where p is an affine parameter of horizontal curves on SU(B a (i)) which project to geodesies on 
B a (i) starting at i. The matrices t A B £ SU(2,C) can be regarded as "angular coordinates". 
Crucial for the applications is that, T° = 7r _1 (i) = {p = 0}, the fibre over i, is such that 
T° ~ SU(2,C). Thus, working on C a (i) provides a representation of B a (i) where i has been 
"blown up" to a sphere. Note, however, that while B a {i) is 3-dimensional, C a is 4-dimensional. 
This extra dimension can be interpreted as giving rise to the spin-weight of the diverse spinorial 
objects when lifted from B a (i) to C a . This issue will arise again in section [71 
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Figure 1: 



Figure 2: The region of the — unphysical spacetime — near spatial and null infinities. The con- 
formal factor leading to this representation is given by = u>fl^ 2 (l — t 2 ). 



The conformal structure of the spacetime (M , g) is exploited by the introduction of conformal 
Gaussian coordinates based on a certain type of conformal invariants: conformal geodesies - 
see e.g. [21]. Accordingly, the spinor field taa 1 will be thought of as being tangent to timelike 
conformal geodesies which are orthogonal to B a (i). For a gauge based on conformal Gaussian 
coordinates, there exists a canonical conformal factor, 9 — i.e. a choice of representative of the 
conformal class of the metric g M „. It is given by 

/ K 2 T 2 \ 2fl 

e = K- 1 n[i , with co = ^ =p + o(p), (4) 

where fl = 0(p 2 ) is the conformal factor of the initial hypersurface — usually obtained from 
solving the constraint equations — , k — n' p, with k' smooth, fc'(i) = 1, and t an affine parameter 
of the conformal geodesies. The function k expresses the remaining conformal freedom in the 
construction. For the ease of future calculations the choice 

K — u>, (5) 

will be made. As it will be seen, with this choice, the locus of null infinity corresponds to r = ±1. 
Other choices are possible, and indeed, the calculations in [24J use the apparently simpler k = p 
which renders a more complicated expression for the locus of null infinity. 

In order to complete the construction of a representation of Af = J + (B a (i)) U J~ (B a (i)) — the 
region of spacetime close to null and spatial infinity — consider the set C QiK = k 1 / 2 ^ of scaled 
spinor dyads {k}/ 2 oa, k 1 / 2 ^}- Define the bundle manifold 

M a , u = {(T,q)\qeC atU ,-l<T < 1}, (6) 

see figure 1. The conformal factor Q given by ((H) can be used to define the following relevant 
sets: 

J= {(r,q) eMa.u I p(g)=0, \t\ < 1}, 

l ± = {{r : q) G M„ | p(q) = 0, r - ±1}, 

= {(r,q) eM a , u | qeC a ^\l°, t = ±1}, 

which will be referred to as, respectively, the cylinder at spatial infinity, the critical sets where 
null infinity touches spatial infinity and future and past null infinity. 

Consistent with the structure of the bundle space M. a .u> it is possible to introduce a calculus 
based on operators X+, X- and X — related to a basis of the Lie algebra su(2,C) — and the 
derivatives d p , d T . For this, use will be made of a frame caa' and the associated spin connection 
coefficients T aa'BC of V. Conforming with the use of space-spinors one writes 

EE' B B' fo \ 

CAA'— TAA'T CEE'—T A'CAB, CAB = T> A C B )B' , (o) 
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with 

T EE 'c EE > = °ab = c AB d T + c\ B d p + c\ B X+ + c AB X_. (9) 

It is also convenient to introduce an unprimed version of Taa'bc, 

^abcd — t b b TaB'CD, Faa'cd —^abcdt b A i, (10) 
which will be decomposed as 

FABCD = —^{£,ABCD - X(AB)Cd) ~ ^AbIcD, (H) 

where xabcd agrees on B a (i) with the second fundamental form of the initial hypersurface. 

The discussion of the curvature will make use of the spinor Qaa'bb'i related to the Ricci 
tensor of the Weyl connection which can be constructed from V and a spinor f aa< with spatial 
part given by f ab — fAA' is related to the acceleration of taa 1 ■ For later use, define 

Qabcd — t c a t d b Qaa'bb 1 - (12) 

The remaining part of the curvature is encoded in the rescaled Weyl spinor 4>abcd — the spinorial 
counterpart of the rescaled Weyl tensor, df 1 vXp = _1 C /1 J/A . 

Proposition 1. There is a frame caa> on M. a .u> for which 

t aa 'c A a> = V2d T , (13a) 
r AA 'f A A'=% fAB = -T CC "T cc , AB , (13b) 

T BB '®AA>BB> =0, e ABC C = 0. (13c) 

The description of AT in terms of the manifold M. a ,un with the conformal factor given by 
and the frame caa 1 of proposition [I] will be called the Friedrich gauge — F-gauge for short. 

In many parts of the sequel it will be necessary to lift objects defined in the base manifold 
J\f to the manifold M. a ,w Scalar and space spinorial objects can be lifted in a natural way, with 
their "angular dependence" expressed in terms of the matrices t A B € SU(2,C). Functions from 
SU(2, C) to C can be expanded in terms of certain functions T n k which can be thought of as the 
lifts of spin- weighted spherical harmonics — see [201 [53] for a more extended discussion: 



In 4- 1 

„ y 4?r T 2n n _ 8 . (14) 

The operators X + , X- and X can be defined by their action on T n l k . The T n fe 's form a complete 
set of orthonormal functions with respect to the Haar measure of SU(2, C) — this fact will be 
used several times. 



4 Conformally flat data without logarithmic singularities 

A discussion of conformally flat initial data which is well suited for the purposes here presented 
has been given in [13] . Its translation into the language of space-spinors has been carried out in 
[53] 154] . Here, the main results required are presented. 

The space-spinorial versions of the momentum and Hamiltonian constraints are given by 

D AB i; A BCD=0, D AB D AB '&=l^ABCDi> ABCD ^- r . (15) 

The spinorial translation of the approach adopted in [T3] — see e.g. [SD] — suggests a solution of 
the form — see the appendix for a definition of the irreducible spinors involved — : 

i>ABCD = -^t(3xabXCD + h ABC D) + %(xABlJCD + X C DUAb) + \{x AB Z C D + XCDVAb) 

pi pi pi 

+ 2^{yABVCD) + 2^(z ab z C d), (16) 



(i 



where 



wit 



Z = X 2 _\* + A + P Q, (17a) 
Tji = -2pd p X-\% + X_\i + pX+Q + iX+J, (17b) 
p 2 = 2pd p (pd p Xi) + X+X-X* - 2Af - 2Af - pdpXi, (17c) 

\ 2 = Xl\ = Xl\ R + XlX 1 = \% + \ ! 2 , (18) 

and A is the lift of a C°° function on B a (i)\{i}. Note that an overbar denotes complex conjugation. 
Finally, 

J = 3(iJx - J 2 )T 2 \ + 3iJ 3 T 2 \ - 3(iJi + J 2 )T 2 \ 1 (19a) 
Q = ^(iQ 2 + Qi)^ ! + ^(Q 3 )T 2 \ + ^(iQ 2 - Qi)T 2 \. (19b) 

In the above expressions, A, denotes a constant related to the expansions conformal Killing vector. 
The term J expresses the angular momentum of the initial data, while Q is associated with the 
so-called conformal boosts. Note that boosted initial data sets — i.e. with linear momentum — are 
excluded from the class under consideration. If the initial data is such that it has a non- vanishing 
linear momentum, then it can be shown that the asymptotic expansion of i? around infinity will 
contain logarithmic terms already at order 0(p 2 ). As firstly noted in [301 [SU], even completely 
regular initial data sets — e.g. time symmetric, conformally flat — may give rise to obstructions 
to the smoothness of null infinity. It is expected that initial data with logarithms will render 
further obstructions. This kind of obstructions arising from particular non-smooth features of 
the data have not yet been studied. Thus, attention is restricted to initial data sets which are as 
regular as possible. 

Even if the initial data sets are non-boosted, there is still the possibility of having solutions 
to the Hamiltonian constraint whose asymptotic expansions contain logarithms. The conditions 
under which the solutions to the Euclidean momentum constraint give rise to logarithm- free 
solutions of the Hamiltonian constraint have been spelled in [13]. From their results it follows 
that if 

A= — + A^, (20) 
P 

where A and are the lifts of C°° functions on B a (i), then the asymptotic expansion of the 
conformal factor i9 contains no logarithms. Recalling that the lift of an analytic function / to C a 
admits an expansion of the form 

oo 2q 

/ = EE Ew^V- ( 21 ) 

P=0 g£Q(p) k=0 

where Q(p) — {0, 2, ... ,p} if pis even and Q(jp) = {1, 3, . . . ,p} if p is odd, the following assumption 
is made: 

Assumption 1. The class of initial data (S, h a p 1 Xafi) under consideration is conformally flat 
and with a solution to the momentum constraint given by H6\) where the lift of the complex 
function A to C a is of the form 

oo p 2q 

a = EEEw^V- (22) 

p=2 q=2 fc=0 

For the class of initial data comprised by assumption [U the lift of the conformal factor d to 
M. a ,ui admits the parametrisation 

i? = 1/p + W, (23) 



Tn older accounts — see e.g. [331 — . \§ and \' 2 are referred to as, respectively, the electric and magnetic parts 
of A... 
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where 

oo p 2p 

W/ = ? + E^^' ^EE"'* 2 ** 3 ^' (24) 
with m the ADM mass, and w p -2 q ,k € C satisfying reality conditions so that T4 7 = W. 
4.1 A further regularity condition 

Attention will be now given to a different type of regularity condition on the initial data. In 
[2"0] it has been shown that, generically, a specific class of obstructions to the smoothness of null 
infinity will be present unless certain coefficients in the expansion of the spinor 4>abcd on the 
initial hypersurface, <S, satisfy a particular relation. In what follows, objects expressible in terms 
of initial data which conrol the regularity of the development at null infinity will be referred to as 
obstructions. The obstructions described in [20] are related to the structure of the principal part 
of the evolution equations (|30b[) , and are shared by zero-rest-mass fields — see gB] . For the case 
of time symmetric data sets these conditions on the Weyl spinor can be rewritten as a condition 
on the Cotton-Bach tensor. Namely: 

d (a p b p ■ ■ ■ D AlBl b AB cD){i) = 0, P = 0, 1,..., (25) 

where bABCD is the spinorial version of the Cotton-Bach tensor. The way the above formula has 
to be modified in order to accommodate initial data sets which are not time symmetric is not 
known. However, the conditions on the coefficients of the expansion of the Weyl tensor given in 
[20] can be calculated explicitly up to any desired order for conformally flat initial data sets. One 
finds that — cfr. with [53j — : 

Ap;2p,fc + \;2p,k = 0, (26) 

with p = 2, 3, • • • , k = 0, . . . , 2p. 

An explicit calculation shows that at least for p = 2, . . . , 6 the conditions (|26p are equivalent 

to 

D {Ap _ lBp ^---D AlBl C% CD) (i) = 0, p=l,...,6, (27) 

where 

Cabcd = V%D E \ a Xbcd)ei (28) 

Xabcd = ^ 2 *I j abcD' with iPabcd being the part of the general solution to the Euclidean mo- 
mentum constraint which is obtained arising from the real part of A, that is, A^ Q. It is noted 
that for the class of data under consideration, C^ BCD (i) — 0, hence this particular condition 
poses no restriction on A. The equivalence of the conditions (|2"6"]) and (|2~7[) to all orders p will be 
discussed elsewhere. 

5 Expansions near spatial infinity 

The evolution of the spinorial quantities c^ B , Fabcd, Qabcd and 4>abcd describing the geome- 
try of the manifold M. a ,u will be given by means of the conformal propagation equations |20U23j . 
The procedure described in section [3J by which the point i is replaced by the cylinder at spatial 
infinity X, renders regular data and equations for the unphysical spacetime. Moreover, it leads to 
an unfolding of the evolution mechanism which allows to analyse the process at arbitrary order 
and in all details. This kind of analysis is made possible by the total characteristic character of 
the cylinder at spatial infinity with regard to the conformal propagation equations — X can be 
understood as a limit set of incoming and outgoing light cones. 

4 The tensorial equivalent of Cabcd ls given by 
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In order to briefly discuss the equations implied by the extended conformal equations it is 
convenient to write: 

v = Kbi Fabcd, Qabcd), (29a) 

= (00, 01,02>3>4), (29b) 

In this manner, one can write the propagation equations concisely in the form 

d T v = Kv + G(v, v) + L<f>, (30a) 
V2Ed T + A AB c%d^ = B(T AB c D )<t>, (30b) 

where K, G denote, respectively, linear and quadratic functions with constant coefficients, L de- 
notes a linear function with coefficients depending on the coordinates via and related objects, 
while E denotes the (5 x 5) unit matrix, A AB c AB are (5 x 5) matrices depending on the coordi- 
nates and B{T abcd) is a linear matrix- valued function of the connection Tabcd with constant 
coefficients acting on 4>. The reader is referred to [20] for a full list of the equations. 

The evolution equations (|30a[) and (|30b[) are to be supplemented with initial data which can be 
readily calculated once an initial data set (5, h a /3, Xap) satisfying the vacuum constraint equations 
is known — see [S2] for the details. 

The unfolding of the evolution process referred to in the previous paragraphs can be thought as 
the calculation of a certain type of asymptotic expansions for the diverse field quantities grouped 
in the vectors v and (j>. Consistent with the class of initial data under consideration — which is 
expandable in terms of powers of p — one can put forward the following expansion Ansatz: 



w ~ E - v( V> ( 31a ) 
~ E ( 31b ) 



The upper limits of the above summations have been intentionally left open to emphasise some 
convergence issues which will be discussed later. Transport equations for the above equations can 
be calculated by differentiating the propagation equations (|30a|> - (|30b|) and then evaluating at X 
— that is, p = 0. These equations are schematically of the form 

d T « W = Kv&> + G(v {0) , ) + G{v^ ,«(°)) 
p-i 

+ E (32a) 
(V2E + A AB (c AB )^)4>^ + A AB (c c AB )^d c ^ = B(Tf BCD )^ 

+ E ( P ) { B ^abcd)^ p ~' ]) - A AB {c%)^d^) , (32b) 
j= l w 

These equations will be referred to as the transport equations of order p. A crucial feature 
of the transport equations is that their principal part is universal. The solution for p = 
is also universal — that is independent of the initial data — and agrees with the solution for 
Minkowski spacetime. For p > 1 they are linear differential equations for the unknowns of 
order p. Furthermore, note that the subsystem (I32aj) consists only of differential equations. The 
transport equations (|32a|) - (|32bp are decoupled in the following sense: knowledge of 0"', 
j = 0, . . . ,p — 1, together with the initial data v^p)\ x o allows us to solve the subsystem (I32a[) to 
obtain the quantities v^ p \ With v^ k \ <^>W, fc = 0, . . . ,p and I — 0, . . . ,p — 1 and the initial data 
<^ p )|xo at hand one could, in principle, solve the equations (|32bj) to get . 

In practise, in order to carry out the above described algorithm one makes use of the fact that 
all the field quantities appearing in the evolution equations (I30a[) - (|30b|) have a definite spin- weight 
and expansion type. Thus, the coefficients r', < j, I < p can be rewritten in terms of the 
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functions T n m k and r-dependent coefficients. For example, for the coefficients appearing in the 
components of the rescaled Weyl tensor one writes: 



— , p 2q 1 

^ = J2 J2 ^-y a 3-.P^( T ) T 2 q k q ~2 +J P P - (33) 
p>|2-i|g=|2-j|fe=0 P ' 

The products of functions T n m k which naturally arise in ()32all - (|32bj) are linearised by means of 
the use of the Clebsch-Gordan coefficients of SU(2, C) — note that the products appearing in the 
principal part of the equations are somehow trivial as they involve the function T ° . The above 
"decomposition in terms of spin-weighted spherical harmonics" allows for an implementation 
in a computer algebra systenf) whereby the diverse transport equations in the hierarchy and 
their solutions are explicitly calculated to a given order p. The order p is restricted by the 
computer resources available — generally up to p — 7 for a spacetime without axial symmetry. 
The assumption of axial symmetry allows to carry the calculations up to a couple of orders more. 

An important question in the programme is the sense — if any — in which the expansions (|31a[) 
and (|31b[) approximate the solutions of the equations (|30ap - (|30b[) . In particular, one would like 
to be able to estimate the residues 



R p+l (v)=v-J2^ iP ' ) P P '> (34a) 
p'=0 p ' 

^ww^-E^V- (34b) 

p'=0 p ' 

The estimation of the above residues and proving that the asymptotic expansions actually do 
approximate the solution to the conformal field equations is one of the remaining outstanding 
problems in the analysis of the structure of spatial infinity. It is expected that some suitable 
extension of the estimates discussed in [22] in the case of the spin-2 massless field (linearised 
gravity) will permit to remove this hurdle. Some of the technical complications of the problem 
lie in the degeneracy of the principal part of the equations (|30bjl and (|32bj) at the critical sets 
Z^, which precludes the use of standard methods of the theory of partial differential equation. 
Further, as discussed in [20l [23l [50l [52 , 53 , this degeneracy translates in a non-smooth behaviour 
of the solutions of the transport equations — which presumably will also extend to the actual 
solutions of the field equations. It can be seen that for the data of assumption [T] and if the 
further regularity conditions ([26]) hold, then up to p — 4 the coefficients in pia|) and (I31b|) have 
all polynomial dependence in r. Starting with p — 5, logarithmic terms which are singular at 
t = ±1 appear. This adds further complexity to the analysis of the asymptotic expansions. 
Hence, most of the calculations of expansions are stopped once the first logarithmic divergences 
have appeared. 

In order to proceed with the analysis of this article, the following assumption will be made: 

Assumption 2. For the class of initial data under consideration, the propagation equations 
L30a)) - l30b}) have a unique solution in M. a ,w Further, the expansions \31a\) and \31 b\) where their 
coefficients are given by the solutions of equations $32a)) and i32b\) approximate these solutions in 
the sense that the residues jj34a\ ) and (fff^fo[ ) can be suitably estimated and are finite in the region 
of interest. The regularity of the solutions is at least that of the r-dependence of the expansions. 



6 The Newman-Penrose gauge 

This section discusses a certain class of gauge conditions on the spacetime (M,g^ u ) near each 
component of null infinity — i.e and — known as the Newman-Penrose gauge — NP 
gauge. This class of gauge conditions is convenient for the analysis of radiative properties of 
the spacetime, like its peeling behaviour and its mass loss. The NP gauge consists of certain 
requirements on the conformal gauge of the spacetime, together with adapted coordinates and 

5 The computer algebra system maple V has been used for this purpose. 
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BC = 00 


BC = 01 


BC = 11 


AA' = 00' 


K 


e 


7T 


AA' = 10' 


P 


a 


A 


AA' = 01' 


a 


a 


A 4 


AA' = IV 


T 


7 


V 



Table 1: The correspondence between the spin coefficients r^A'sc and the standard NP notation. 

frame. The NP gauge conditions have to be implemented independently for and . Here, 
for conciseness the discussion is concentrated on / + . The analogous construction for can 
be obtained, in principle, in the same way by interchanging the role of the relevant vectors. The 
construction of the NP gauge is closely related to the so-called asymptotic characteristic initial 
value problem in which initial data is prescribed on an outgoing null cone and a portion of 
and allows to recover a portion of the spacetime to the past of these null hypersurfaces — see 
[T71 IT51 [TBI, [T9l 129"] . Alternatively, one could prescribe data on an incoming null cone and on a 
portion of J?~ . Here, however it will be assumed that the spacetime has been already constructed, 
and a suitable conformal factor, 0, has been found, allowing to locate the conformal boundary 
of the spacetime. 

Let g^v = Q 2 g^ v denote the metric of the conformally rescaled (unphysical) spacetime. Our 
discussion will consider frame fields (e-AA' Y such that g(eAA' , cbb') = (ab^A'B 1 - The correspon- 
dence to the standard NP notation is given by (eoo') M = l 1 *, {ewY = n M , {eoi'Y = mfJ " and 
(eio'Y = '■ The spin connection coefficients associated to the frame (caa'Y are given by 

Taa'BC = \{{eAA'Y{eBvYV„{ecw)v + {e A A>Y{e cV yV ^e B0 ,) v ). (35) 

The correspondence with the standard NP notation can is shown in the table 1. 

Consider a subset M + C M. containing a portion of J^ + . And let ^ + be a fiduciary cut of 
null infinity: <*f+ C J? + , c € + C J\f+. It shall be assumed that J?+ U N + ~ Sf+ x (-e,e), e > 
— that is, it contains an open set of cuts. Note that because of the assumption that the spacetime 
M. arises as the development of some asymptotically Euclidean initial data set, it follows that 
the cuts of have the topology of S 2 . 

The following result shows that there is a certain gauge adapted to the description of J'^ . 

Proposition 2. Assume thatj\f + is suitably small. Then, there exists onAf + a conformal factor 
0+, a frame {c\ A 'Y an( ^ coordinates (tt+, r+, x 2 ,, x+) such that 

LlO'OO = ^Ol'O'O'' Ifi'oo = ^Ol'Q'l' + ^01'01> ^OO'AB = 0> on (36) 

and 

^io'ix = -^li'ii = -^oi'ii = r xl , 01 = r 11/00 = r 10 , 00 = o onj' r\M , (37) 

where ^AA'BC are ^ e connection coefficients with respect to the frame {c\a'Y °f ^ e Levi-Civita 
connection V + of the metric g~^ v = {Q + ) 2 g^u, where + = 6 + Q. Furthermore, 

R+ =0 on J\f + (38) 

and 

= \ Utu \C R+ ) (4vY(4vY = 0on,y+ nM + , (39) 

where R + and Rz v denote, respectively, the Ricci scalar and Ricci tensor of the metric g~^ v . The 
coordinate u + on J? + is an affine parameter of the generators of null infinity. The vector (e^,)^ 
on is tangent to the generators of null infinity and satisfies e^/(ti+) = 1. OnM + , u + satisfies 
the eikonal equation V +M u-|-V+'M+ = 0, and the vector {&qqiY is tangent to the generators of the 
hypersurfaces IA U+ = {u + = constant] and r + an affine parameter. The coordinates x 2 , and x\ 
are chosen such that the metric induced on the cuts of is the standard S 2 metric. 
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Remarks. 

1. Proofs of the above statements can be found in the discussions of [24l [29] . 

2. The conformal factor, frame and coordinates given by the above result will be referred 
collectively to as the Newman- Penrose (NP) gauge for ^ + . 

3. A corresponding NP gauge consisting of a conformal factor 0_, a frame (e AA ,)^ and coordi- 
nates (v- , r_ , x 2 _ , x^L ) can be obtained by performing the interchange i— > 1 , 1 i— ► in the 
spinorial indices, so that e.g. rg 1/Q0 = on CiJ\f~ for some suitably small Af~ C M. 
Here T AA , BC denote the spin coefficients with respect to {(^aa')^} °f the Levi-Civita con- 
nection V~ of the metric g~ = (Q~) 2 gfj, v , with 0~ = 0_0. 

4. The only non- vanishing coefficients on J? + C\J\f + are r^gg, rg" 1/01 and r+ , 01 . In particular 
c + = r^Q'oo cnc °dcs information of physical interest. 

7 Relating the F-gauge and the NP-gauge 

Having introduced two different gauges for the analysis of the class of spacetimes under consid- 
eration, it is now necessary to discuss the way they are related to each other. The NP-gauge is 
hinged on null infinity and is adapted to a characteristic initial value problem, while the F-gauge 
is based on a Cauchy initial value problem. In the F-gauge the conformal factor O, and hence the 
location of null infinity is known a priori by inspection of the initial data on S — as discussed in 
section ([3]). Further, because of the asymptotic expansions of section [5] which allow to examine 
the evolution of the field up to a given order p in terms of Cauchy initial data, the setting on 
which the F-gauge is implemented will be regarded as the starting point. 

As in the previous section, in order to fix ideas, the discussion will be mainly concerned with 
the NP-gauge associated to J? + . A first important observation is that the construction of the F- 
gauge and the associated regular Cauchy initial value problem near spatial infinity — see [20l [50l 
I52[l53j — lives in the 5-dimensional bundle manifold M. a ,ui, and not in the spacetime manifold, M.. 
The correspondence between objects on M. a ,u and M. is carried out by the projection 7r. Clearly, 
the correspondence between objects in the two manifolds is not one to one. This ambiguity is 
related to the spin-weight of the diverse spinorial objects under discussion — put in other way, 
when working on A4 a ,u> one works with all possible choices of frame vectors (eoi')'*, (eio') M - 

Assume that the region Af + considered in section [6] is such that Af + C Tr(M a ,uj)- In order to 
relate the frames caa' — living on M. a ,u — and e A A' — living on M + — it is convenient to fix a 
section of M. a ,u> — to be thought of as a choice of a tetrad c* AA , on Af + . Let U be an open subset 
of S 2 . Let s :U -> SU{2) be a local section of the Hopf fibration SU{2) -> SU(2)/U(V). The 
map s induces a smooth section U x R x E 3 (p,r,p) (s(p),T,p) £ M a ,iu- The image of s 
will be denoted by w . The vector fields tangent to s(U) which have a projection identical to 
that of X± are of the form X± + a±X with some functions a± on s(U) which satisfy a_ = — a + . 
Due to the commutation relations obeyed by X± and X , the functions a± do not vanish on open 
subsets of s(U). Thus, the fields c* AA , on M* a u which satisfy tt(c* aa ,) = n^(c AA ') are given by 

C *AA' = c aa' + (a+c% A , + a-c AA ,)X. (40) 
The spin connection coefficients associated to the frame c* AA , defined on M* a w are given by 

F*aa>bc = ^aa'bc + {a+c\ A , + a-c AA ,)(e OB eco - eiseci)- (41) 

In what follows, objects defined on Af + will be lifted to the section A4* K . In order to avoid 
complicating the notation unnecessarily, lifted objects will be designed with the same symbols as 
their original version on M + . From the context it should be clear where the diverse objects live. 
Let 

n A 'BC = ^ (r + tA(BC* c)A > (In 6+)) , (42) 

the spin connection coefficients associated to the frame c AA , = O^c*^,. From the discussion in 
|2"4] it follows: 



12 



Proposition 3. There is a Lorentz transformation, A A B £ SL(2, C) such that 

e\ A , = A B A A B ' A ,c BB , = -±-A B A A B ' A ,c BB , , (43) 

where e\ A , and 9 + are the frame and conformal factor of proposition [H 
Remarks. 

1. The matrix elements A A B on J? + can be determined entirely by means of a construction 
intrinsic to null infinity. Furthermore, for the class of spacetimes under consideration 

A°iU+ = p 1/2 + dV /2 ), A\U+ = \x+Wi +<IV /2 ), (44a) 
A°oU + = -^X-Wip^ + 0{pV% A\U + = -p-V* + Dip 1 ' 2 ), (44b) 

where W\ is the O(p) term in the expansion of W — see equation (|24|) . Latter discussion 
of the Newman-Penrose constants will require the knowledge of the quantities eQ ,(A A B ). 
One has that 

4a'(A B c) = -A F A A F ' A ,A H C T* FF B H + A B D T\ A , D C . (45) 

Remarkably, the gauge conditions (f37|) enable a calculation of e^ , (A A B ) completely intrinsic 
toJ^+. 

2. The determination of the values of 0+ and d8+ on J? + — which is equivalent to the knowledge 
of the quantities ej^, (#) — can also performed in a way which is intrinsic to 

3. Consistently with the discussion of the previous section there is a function, /+, on such 
that 

(e+,) M = /+V +M 9 + on J?+. (46) 

4. The retarded time it+ on has the expansion 

u+\s+ =\/2^-^ + 4mlnp + u*+ (^rn 2 - 3A + y wA p + 0(p 2 )^ , (47) 

where it* is an arbitrary function of the angular coordinates associated to the supertrans- 
lational freedom. This expansion has been calculated up to 0(p 4 ) inclusive. The retarded 
time V- has a similar expansion on 

v-\s- =V2Q-4mm ( o + u*- ^m 2 +3i + y» r 1 j/) + 0(p 2 )j . (48) 

In the sequel — and in particular for the evaluation of the NP constants — it will be natural 
to push the fiduciary cuts of null infinity to the critical sets where null infinity touches spatial 
infinity. Thus, the following assumption is made. 

Assumption 3. The fiduciary cuts and of, respectively, and J*~ can be made 
coincide with the critical sets 2 + and T~ . 



8 Peeling properties 

In view of the existence of obstructions to the smoothness of null infinity of the developments of 
the class of initial data sets of assumption Q] discussed in sections (|4.1[) and ([5]) it is natural to ask 
how do they affect the asymptotic behaviour near null infinity. The translation between gauges 
discussed in section [7] provides an adequate way of performing this discussion. 

The Peeling Behaviour is a specific type of decay of the components of the Weyl tensor with 
respect to a frame satisfying the NP-gauge. The classical discussion of the peeling behaviour 
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is generally carried out in the physical spacetime — a manifold without boundary. Let g^ v the 
physical metric obtained from g+ via 

9»u = f 2 g+ v , (49) 

where f = l/r + . Coordinates in (M^g^) can be introduced by considering (f, 

Recall that r + = — i.e. f — > oo — gives the locus of J? + . A frame, baa 1 , and a spinor dyad, 

{o+, f4} on the physical spacetime can be introduced via 

e^=^'- 2 e+ A „ (50a) 
5^ = f^oi, 1% = 4. (50b) 

In the coordinates (f, u+, 2;+, x+) one has that Bq , = df. 



The original Peeling Behaviour described in the classical literature assumes sufficient regularity 
of the solutions of the Einstein field equations at, say, future null infinity so that in the physical 
spacetime the components of the Weyl tensor admit an expansion of the form — see e.g. [31 HOI 
[33 [38]: 

* = O(r 5 ), *i=0(f 4 ), ^ 2 =0(r 3 ), V 3 = 0(r 2 ), $4 = O(f). (51) 

where ^0 = ^ abcdo+o+o+o+ , ^1 = ^ abcdi+o+o+o+ , etc. The decay conditions for the 
components ^f n are not independent. Assuming a particular decay for "Jo, it is possible to deduce 
the decay of the other components. Note that smoothness at null infinity is a much stronger 
requirement than that of peeling behaviour for which only a finite degree of differentiability is 
needed. 

In order to proceed some observations of generic nature are made. It is recalled that the 
rescaled Weyl spinor, 4>abcd, is related to the (conformally invariant) Weyl spinor ^ abcd by 

^ABCD = Q&ABCD, (52) 

so that if ^abcd = 0(@), then §abcd = O(Q ). Further, the components ^>abcd are related 
to those in the physical manifold via: 

*„ = e 4 -"*„. (53) 

Moreover, it is noted that 

= O(l/f), (54) 

close to J? + . 

Using the above relations, and assuming that assumption [2] holds, it is possible to read off 
more or less what the Peeling Behaviour of the developments of the initial data sets of assumption 
[T] should be. It is nevertheless necessary to analyse whether there occur any cancellations which 
may improve the decay of the Weyl tensor. 



8.1 Data with D (Ap _ li3p _ 1 • ■ ■ D AlBl C ABC D)(i) = for p = 1,2,3,4. 

The best behaved case that will be considered here is that of an initial data set where the regularity 
condition (|27j) holds up to p — 4 inclusive. In this case, the main result of [52] shows that all the 
functions aj tP -2 q ,k(T) appearing in the expansions (|3"3"]) up to p = 5 are polynomials in r save for 
the functions Oj 5 ; 4 fc (t) , j = 0, 1, 2, 3, 4, with k = 0, . . . , 4 which contain logarithmic terms. More 
precisely, one has: 

a ilB ;4,k = T 5 ;4,fc ( (1 - t-) 7 - j Pj(t) ln(l - r) + (1 + t) 3+ ^P 4 ^(t) ln(l + T )j + Q(r), (55) 

where T^.^k are certain coefficients (obstructions) completely specifiable in terms of the initial 
data, Pj{r) denotes a certain polynomial of degree j in r such that Pj(±l) ^ 0, and Q(r) is an 
arbitrary polynomial. 
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Now, it is noted that from formula (j4|) for the conformal factor 9 it follows that close to, say 

^ = (TT^ a (56) 

Further, (1 + r) _1 is analytic close to r = 1. Hence, a direct use of the relations ([5"5)l . (|52p and 
(|53[). shows that for this class of developments: 

*« = (j^r) • ( 5? ) 

for n = 0, . . . , 4. The first non-smooth terms in the expansions of ^> n arise in \to at order C(l/f 8 ), 
where a term of the form 

T 5 ;4,fcr- 8 lnf, (58) 

k = 0, . . . , 4, is to be found. It is mentioned by passing that in this case the obstructions Ts-^k 
correspond to the polyhomogeneous constants of motion discussed in [5J 031 [Hj . 

8.2 Data with D (Ap _ lBji _ 1 • ■ ■ D AlBl C ABCD )(i) = for p < 3. 

As discussed in [SOI [S5], the obstructions Ts^fc are special in the sense that they are the first 
truly non- linear obstructions — that is, they arise as a product of the specific non-linear features 
of the conformal Einstein equations. In contrast, the regularity conditions (|27p arise from features 
which are present already at the linear level. Nevertheless, it is of interest to analyse how these 
kind of obstructions affect the peeling behaviour of the components of the Weyl tensor. 

As discussed in if X P] 2 P ,k + \-,2 P ,k = 0, for a certain p, then the functions aj jP; 2 Pl fe(r) in 
the expansion (|33[) contain logarithmic terms. More precisely, 

a M ,k(r) = A P ;2 P ,kUl-rr- 2+j (l+r) p + 2 ^' ln(l-T) + (l+r)P- 2 +^(l-r)P+ 2 -J ln(l+r)^) +g(r), 

(59) 

with the obstruction A p; 2 Pl fc such that A p -2 Pt k = if and only if X p -2 P ,k + X p -2 P .k = 0. As discussed 
in section |4~T1 for p = 1, . . . , 6 one can explicitly show that it is equivalent to the condition ([2"T)l . 
It is plausible that this kind of logarithms will propagate — due to the hyperbolic nature of the 
propagation equations — and may produce further non-smoothness at higher order. This process 
is, however, not yet understood. 

According to the above discussion, if for p — 4 the regularity condition ([27| does not hold, 
then the resulting peeling spacctime will contain — in addition to the logarithmic terms of ([58"f 
the following singular terms: 

A 5; io,fc^ 8 hif. (60) 

Note that in contrast to Ts^^, which is an object of quadrupolar nature, A^-iq^ is of 2 5 -polar 
nature. 

The above argumentation can be easily extended to the case where the obstructions A p; 2 Pj fe 
with p = 2, 3 are present. In these cases, the component of the physical Weyl tensor contains 
singular terms of the form 

A p , 2p .kf- 3 - p lnf. (61) 

It can be checked explicitly that when performing the Lorentz transformation implied by the 
matrix A A B £ SL(2, C) of section [3 there is no cancellation taking place which would improve 
the regularity of the components of the Weyl tensor. Thus, in the case p = 2 the singular terms 
arise sufficiently low in the expansions to preclude the peeling behaviour. 
One obtains the following. 

Proposition 4. Under assumption^ the developments of initial data sets of assumption^ with 

D {AlBl C ABC D)(i)=0 (62) 

in the sense that 

* = C(r 5 ), *i=C(f 4 ), * 2 =C(f 3 ), * 3 = C(r 2 ), * 4 = e>(r). 
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On the other hand, if 

D { a iBi Cabcd){*)^0 (63) 

then 

*o = C(f~ 5 lnf), (64) 
and accordingly the development does not peel. 
Remarks. 

1. Note that the requirements on conformally flat initial data to have peeling are very mild. In 
particular initial data sets like the Brill-Lindquist 5 , Misner [32] and Bowen-York [4] will 
satisfy them. Thus, modulo assumption [2l the developments of these initial data sets will 
peel. However, the Weyl tensor of the unphysical spacetime will not be smooth due to the 
presence of the terms in (f58| — the obstruction Ts^^ has been explicitly been calculated 
for this case in [T3] and shown to be non-zero unless the data is Schwarzschildean. These 
properties have been used in [51) to show the non-existence of conformally flat slices in the 
Kerr spacetime. 

2. On the other hand, as shown in [llj . the Einstein field equations are compatible with 
characteristic initial data with a decay of up to = 0(lnf/f 3 ). The question is, which 
kind of Cauchy data would produce such type of behaviour? From the discussion in this 
work, it is to be expected that conformally flat initial data with linear momentum and/or 
for which the complex function A has a singular behaviour worse than that of (|20[) will give 
an example of initial data which would render this kind of spacetimes. 

3. The definition of the Newman-Penrose constants — see section [TD] — requires certain degree 
of smoothness of the Weyl tensor. From the discussion presented in this section, the required 
condition on the initial data is 

D (A 3 B 3 D A2 B 2 D AlBl C A BCD){i) = 0, (65) 

so that 

*o = 0{f- 7 \nf). (66) 

For spacetimes whose initial data do not satisfy the above it is not hard to see that the 
obstructions A2 ; 4 i ? £ and A3 ; 6,fc are related to the conserved constants and logarithmic NP 
constants of P [431144], 



9 Expansion of a on null infinity 

As another application, the expansions of the Newman-Penrose coefficient cr + = F^ , QQ on 

— the asymptotic shear — are discussed. The relevance of this stems from three — in principle 

unrelated — directions. Firstly, as discussed in [34], if (a + ) R — the electric part of a + — satisfies 

lim (<t + ) r \jt+ = 0, (67) 

U+ — > — OG 

that is, the limit as p — > 0, then there is a canonical way of selecting the Poincare group out of 
the asymptotic symmetry group — the Bondi-Metzner-Sachs (BMS) group. Secondly, a + \j;+ has 
been called in older accounts the news function, and enters in the expression of the mass loss 
formula. Here and in the sequel, the overdot denotes differentiation with respect to u+. Thirdly, 
as discussed in [6[ 1 1 5 j . cr + |,y+ is related to the memory of non-linear tail waves first discussed by 
Christodoulou [HH]. 

Direct computer algebra calculations with the asymptotic expansions of section [5] and the 
transformation between the F-gauge and the NP-gauge render the following result. 

Theorem 1. For the class of initial data of assumption^ and under the existence assumptions 
[H and [3] one has that 

a+U + =0(p 3 ), (68) 

as p — > 0. 
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Remarks. 



1. On the light of this result it is natural to ask what kind of initial data is required to obtain 
a development such that the decay condition 

cr + U+ as 0, (69) 

occurs. Preliminary calculations suggest that boosted conformally flat initial data would 
give rise to developments for which the decay (|69[) occurs. Another alternative would be 
to move away from the choice (|20j) of the complex function A giving rise to the higher 
multipolcs of the second fundamental form. In any case, it seems that in order to obtain 
the asymptotic decay (|69j) one has to consider solutions to the momentum constraint with 
logarithms in their asymptotic expansions. These issues will be explored elsewhere. 

2. Another natural question to ask is: what is to be expected if one moves away from the 
restriction to conformally flat data? It is to be presumed that as long as one restricts to 
a class of "reasonably smooth" conformal metrics, h a p, the conclusion of theorem [T] will 
still hold. As an example of what one would call "reasonably smooth" consider metrics 
with the same smoothness as the conformal metric of stationary solutions that is, h a p G 
C 2 ' a (B a (i)), where C k ' a denotes the standard Holder space, see e.g. [12] . 

3. The above two observations and the results of section [5] give credence to the conjecture that 
any peeling spacetime will have o~ + \^+ — ► as u+ — * — oo. 

4. If the initial data set satisfies the regularity condition (f2"T|) to some orders, then it is possible 
to provide further information about the expansions of a + on null infinity. In particular, 
for a spacetime such that the first obstructions to null infinity appear at order p = 4 in the 
components of the Weyl tensor, one has that 

= -\{^ X - W ^ - ( X - W i f + |t(A_ J) 2 - Z1XJW\X-J 



+ ^V6m(A 2 - A 2 ) - 4V6(A 3 - X 3 ^jp 3 + 0(p A ). (70) 



From the expansion ([70]) it is possible to derive an expansion for the Bondi mass of the 
development, mg, close to the critical set I + . In the NP-gauge one has that 

m B = -^£ ^0ABCDA A 1 A s 1 A c o A D o)U++^ + l.^+^ + U+)d^. (71) 

The above expression can be lifted in a natural way to the manifold M a ,u>, by noting that surface 
element dS has to be replaced with dS'da where a is a parameter along the fibres. The functions 
T m k l form an orthonormal basis with respect to the Haar measure on SU (2, C) which is given by 
d/i = (l/47r 2 )dSda, thus only terms with T Q ° contribute to the integration. One finds that: 

Proposition 5. For the class of initial data of assumption^ and under the existence assumptions 
[H and [3] one has that 

lim tub — rn. (72) 

That is, the Bondi mass agrees with the ADM mass on the critical set I + . 
Remarks. 

1. Versions of the above result with different regularity assumptions can be found in the 
literature — see e.g. [1][28]. The present form has been anticipated in [24] , 
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2. Using the expansions (|70|1 — which are valid for developments for which the first obstructions 
to the smoothness of null infinity appear at order p = 4 — one can obtain a generalisation 
of the expansion for the Bondi mass given in [47) : 

37r 2 4 
ran = m > \Hu 



where 



fe=0 



3^ 
560 



El** 



1 



h0(l/« 8 ), 



(73a) 
(73b) 



fc=0 



— 15 

-ffo = (nm>2;4,0 ~ 2 ^6 w i-2.o) 7r^>{-h + iJi) 



i?l = (mu>2:4,l 



i?2 = (™W2:4,2 



if 3 = (rnw 2 ;4,3 



i?4 = (TOW2:4,4 



+16V6wi ; 2,o(«/2 - iJi) 

+16m(A 2; 4,o - A 2;4 ,o) - 24(A 3;4 , - A 3;4 ,o), (74a) 
4v / 3wi ; 2 : owi;2a) — 15v / 3( Jl J3 + iJ 2 J 3 ) 
-16V3(iJ 3 wi ;2: o - (J2 - iJ\)wi;2,x) 

+16m(A 2;4 ,i - Aa;4,i) - 24 ( A 3;4,i - A 3;4 ,i), (74b) 

4 ™1;2,1 " 4Wl;2,0Wl;2,2) + 15(jf - j\ - J%) 

-16(iJl(lt)l;2,2 - U>1;2,0) - ^2(^1:2,0 + »1;2,2) + 2iJ 3 «;X;2,l) 

+16m(A 2;4i2 - A 2;4 , 2 ) - 24(A 3;4 , 2 - A 3;4i2 ), (74c) 
4v / 3i0i ; 2,i' ly i;2,2) + 15vo(Ji J 3 - iJ 2 J 3) + 
-16V^(U 3 Wi- 2< 2 - {-h + iJiVi;2,i) 

+16m(A 2; 4,3 - A 2;4 , 3 ) - 24(A 3;4 , 3 - A 3;4 , 3 ), (74d) 
15 



+16V6wi ;2i2 (J 2 + iJi) 

+16m(A 2;4 ,4 - A 2;4 ,4) - 24(A 3;4 , 4 - A 3;4 ,4)- 



(74e) 



In [15] , Frauendiener has identified the non-linear memory effect of tail waves in terms of the 
difference, Aa + \y+, between the value of <J + \.^+ at a cut of J< + given by c to u+ = {u+ = constant} 
and its value at some fiduciary cut c € + . In the current context it is natural to let the critical set 
I + to be this cut. In this light theorem [1] states that there is no "residual" distortion contained 
in the initial data, and that, in principle there is no restriction to the possible values of Aa + \jr+ 
which can be attained along J? + . On the other hand, if a + \^+ -/* as u+ — ► 0, then there are 
values of the memory effect which are not accessible for Aa + \jf+, as it satisfies the equation 



X 2 _A^ 



r+<7 



r+cm 



(75) 



10 The Newman- Penrose constants 

The Newman- Penrose constants are a set of nontrivial complex quantities of quadrupolar nature 
defined on cuts of ^ + (five of them) and J"~ (another five) which are absolutely conserved — see 
[531 [33] . Here, absolutely conserved means that their value is independent of the choice of cut, and 
that they are supertranslation invariant. Further, the two sets of quantities transform adequately 
under Lorentz transformations — i.e. they form a representation of the Lorentz group. 

The argumentation leading to the NP constants makes certain assumptions on the regularity 
of the spacetime at null infinity. Let = ^ abc DO+o+d+d 1 ^ . 

Proposition 6 (Newman-Penrose, 1965). Assume that on J\[ + \ (JV + n the component ^0 
of the Weyl tensor has the decay 

*o = "f^f- 5 +^ ] f- 6 +0{f- 7 \nf), (76) 
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then the quantities (the NP constants for ^ + ) 



G+ = A 2 Y 2m ¥^dS, (77) 

with m — —2, ... ,2 are supertranslation invariant and absolutely conserved — that is, G+ = 0. 
Remarks. 

1. Note that the spacetime is required to have a bit more regularity at null infinity than what 
is needed for the peeling behaviour. 

2. An analogous result holds for if on Af~\(Af~r\t#~) the component ^4 = ^ abcd^+^+1+^+ 
has the decay 

§ 4 = -qi ( pf' 5 + 4<f V 6 + 0{f- J mf), (78) 
where here f = l/r_. Then the NP constants for are given by 

G- = <f 2 Y 2m ¥^dS. (79) 
Jv- 

3. Note that in principle there is no reason to expect a relation between and G~ , unless 
the spacetime is time reflexion symmetric — like in the case analysed in [24] . 

A clear cut physical interpretation of the NP constants has remained, so far, elusive. From 
their evaluation in particular examples and from comparison with linear theories, it follows that 
the content of the constants is nontrivial [TJ] HH [35] [31] . There has been some work attempting 
to relate the constants with the radiation content of the spacetime at late times and the value of 
the Weyl tensor at timelike infinity, i + — see e.g. [35] EH [45] . Remarkably, it has been recently 
shown by Wu & Shang [56] that the NP constants vanish for algebraically special spacetimes. 
This was a long-standing conjecture — see [3U]- From here it follows that the constants vanish 
for the Schwarzschild spacetime (which was known) and for Kerr — i.e. for the (presumably!) 
asymptotic states of spacetimes containing black holes. This result, together with conjectures 
given in [52[ 153] support the idea that the NP constants measure in some way the gravitational 
radiation content of the spacetime at earlier and later times. 

In [21] it has been shown that it is possible to reexpress the NP constants in terms of quantities 
defined on the unphysical manifold, Ad. For this, note that (f>o — \P abcdo+o+o+o+ , so that one 
obtains 

G+ = -^-/ 2 y 2n e+,(^ )d5. (80) 
In order to perform the obvious lift to the bundle space Ad* auJ recall the correspondence 

2 Y 2n ^(-if- n ^T i 2 ' n Q , _ 2 Y 2n ^(-i)^^ Ti 2 -\ (81) 
If d/i denotes the Haar measure on SU(2,C), then G+ lifts to 



G+ = ^ f T 4 2 -V + ,(<MdMda, 



(82) 



where a denotes a parameter on the fibres of the fibre bundle Ad a ,w The second integration 
does not alter the result as the integrand is independent of a. For the NP constants on J?~ one 
obtains in a similar way that 
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Expanding the integrands one finds that 
G+ = ( -^^ ^ £ + T 4 2 "" 4 1 ( A B A c A^ A B (A A A A ' , c^, (0 BCD£; ) - 30 BODB e+ , (0)) 

+4dA B Q A c AVoo< (A £ )^ BCDB ^ d/ida. (84) 

(-i) 2 "" / 5 / „ 2-n 1 



G- = -^-y^f ^-^(AV lAV^M* v c* aa ,^bcde) - WBCDEelvW)) 

+48 A B , h c 1 A^u' (A^J^socbJ dfida. (85) 

The strategy is now to calculate 

G+ | I+ = lim G+ G" | x - = lim G~ , (86) 

in order to obtain an expression for the NP constants in terms of initial data quantities. A 
direct inspection using the results of the asymptotic expansions discussed in section [5] reveals 
that there are potentially conflictive terms of order 0(p~ 1 ). Remarkably enough these cancel 
out. A computer algebra calculation renders the following result. 

Theorem 2. For the class of initial data given by flp and under assumption^ the Newman- 
Penrose constants at are equal to those at . 

More precisely, one has that 

•2-n 

G„| x+ = G n \ T - = —VT^G' n , (87) 

where 

113 

G' = l27(mw 2 . A , - 2V&w\. 2fl ) - — ( J 2 + iJi) 2 
+240wi ;2 ,o(J 2 -iJi) 

+ 120m(A 2;4 ,o - A 2;4 ,o) + 12(A 3;4 , - A 3;4:0 ), (88a) 
113 

G[ = 127(tow 2:4 ,i - 4V3wi ;2 ,oWi ;2 .i) - — vT2(Ji J 3 + iJ 2 J 3 ) 
-120\/2(iJ 3 wi ;2 ,o - (J 2 - iJi)wi ;2 ,i) 

+120m(A 2 ; 4 ,i - A 2i 4,i) + 12(A 3; 4,i - A 3;4 ,i), (88b) 
G' 2 = \27{mw 2 . Aa - Aw\. 2A ) + 113(J 3 2 - J 2 - j|) 

-240(iJi(wi ;2 , 2 - lOl ;2 ,o) ~ ^2(wi ;2j0 + tOl;2,2) + 2iJ 3 l0l;2,l) 

+ 120to(A 2;4j2 - A 2;4 , 2 ) + 12(A 3;4 , 2 - A 3;4 , 2 ), (88c) 
113 

G 3 = 127(touj 2;4 , 3 - 4V3u>i;2,iWi.i2,2) + — yT2(Ji J 3 - iJ 2 J 3 ) + 
-120V2(iJ 3 w 1;2 , 2 - (J 2 + iJi)wi ;2 ,i)) 

+ 120m(A 2;4:3 - A 2;4:3 ) + 12(A 3;4 , 3 - A 3;4 , 3 ), (88d) 
113 

G 4 = 127(mw 2;4 , 4 - 2V6wj. 2a ) + — (J 2 - iJi) 2 
+240wi ;2:2 (J 2 + iJi) 

+ 120to(A 2;4j4 - A 2;4 , 4 ) + 12(A 3;4 , 4 - A 3;M ). (88e) 
It is noted that the integrand of the formula (|84[) can be written concisely as: 

G = mW 2 + V6(X_T^i) 2 + ^-%/6(A_ J) 2 - ^-V&X-W X X-J +^-m{\ 2 - A 2 ) 

4572 127 127 

+ ^(A 3 -A 3 ). (89) 
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From where it can more easily be seen that the constants have the structure: 

(mass) x (quadrupole) + (dipole) 2 + (angular momentum) 2 — (dipole) x (angular momentum) 

+ (angular momentum quadrupole). 

In the view of the above formulae, it is natural to ponder what would happen if one were to 
consider a more general class of initial data sets. It is to be expected that the theorem [2] will still 
hold if one considers non-conformally flat initial data, for the time asymmetry in the solutions of 
the Einstein constraint equations is fed by the second fundamental form, not by the conformal 
metric. If one were to consider boosted initial data, it is quite plausible that null infinity may not 
be regular enough — recall the existence of logarithms in the initial data discussed in section 0] — 
to ensure the constancy of and G~ — unless, perhaps if the initial data set is not maximal, 
like in the case of boosted slices of Schwarzschild. 

11 Concluding remarks 

The present article has been concerned with an analysis of the physical implications of certain 
asymptotic expansions for asymptotically flat spacetimes which allow to relate the behaviour 
near null and spatial infinity to properties of the initial data. These asymptotic expansions are 
obtained in a gauge which exploits to the maximum the conformal structure of the spacetime. 
Hence, in order to gain intuition on the physics of these expansions one has to transform to a 
more "physical gauge" . Crucial feature of the expansions is that they bring to the foreground the 
fact that the development of generic initial data will render a non-smooth conformal boundary 
of the spacetime. This non-smoothness manifests itself in the appearance of logarithmic terms in 
the solutions to the transport equations which generate the terms in the asymptotic expansions. 
It has to be emphasised that these logarithmic terms are conformal invariants of the spacetime 
and not the product of some bad choice of gauge. In other words one could completely forget 
about the F-gauge and the cylinder spacetime, but one would still have to confront in one way or 
another their implications — as the discussion about the peeling behaviour in section [5] evidences. 
On a more practical side, it is not clear what the implications for numerical simulations of the 
non-smoothness of null infinity are. As seen here, as long as no linear momentum is included in 
the reference asymptotic end, the smoothness of most conformally flat initial data sets of interest 
is enough peeling so that most of the formalism remains untouched. 

An interesting and in principle not evident result from our calculations is the fact that the 
NP constants — in the cases for which they are well defined — offer a set of absolute conserved 
quantities for the whole spacetime, and not too different for each component of null infinity. It 
has been argued that this feature will be preserved in the case of more general initial data sets. 
Further, a simple computation also reveals that the logarithmic NP constants should also share 
this feature. 

Finally, it has been seen that under certain existence assumptions having the peeling behaviour 
automatically implies the right kind of behaviour of the asymptotic shear at null infinity to obtain 
a canonical representative of the Poincare group out of the BMS group. This result exemplifies 
the power and insights that can be gained by methods presented in this work. At this point it 
would be natural to proceed with a discussion of angular momentum. However, as the result of the 
non-existence of conformally flat slices in the Kerr and other (non-static) stationary spacetimes 
exposes, there is a certain non-cooperation between conformal flatness and the presence of angular 
momentum at null infinity — see e.g. [5TJ 08]. So it seems that a more natural setting for this 
analysis is that of non-conformally flat metrics. 

This last issue and other questions raised in this article will be further dwell upon elsewhere. 
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A Definitions of some spinors 

In the text, use of the following spinors has been made: 

xab = ~^2>° aLb + l aOb), Dab = ~~/^ LaLb > Zab = ~^° A ° B ' '■ 90a -' 

(°ABCD = °AOBOcOD, CaBCD = °(AObOcI>D), ^ABCD = ^A^B^C^D^ (90b) 

h-ABCD = -£A(C e D)B- (90c) 
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